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Abstract 

A recent scheme for perfect transmission of quantum states through quasi-one dimensional 
chains requires application of global control at regular intervals of time. We study the effect of 
stochastic noise in this control and find that the scheme is robust for reasonable values of disorder. 
Both un-correlated and correlated noise in the external control are studied and it is remarkably 
found that the efficiency of the protocol is much higher in presence of correlated noise. 



1 Introduction 

Since the early developments in the theory of quantum information, the task of coherently transferring 
quantum states through long and short distance communication lines has been of great importance. 
While photons are the ideal carriers of quantum information over long distances, it has become evident 
that the best possible method for transferring quantum information over short distances is to exploit 
the dynamics of many body systems, specifically of regular arrays of qubits constituting suitable spin 
chains. In this framework, an arbitrary qubit state is first coupled to the array and then carried to des- 
tination by the natural dynamics of the whole system, where it can be extracted with certain fidelity. 
This idea was first introduced in |lj, where it was shown that the natural dynamics of a Heisenberg 
ferromagnetic spin chain can achieve high-fidelity transfer of qubits over distances as long as 80 lat- 
tice units. In contrast to this traditional "passive" protocol, different approaches soon emerged. One 
idea was to engineer the couplings between the various spins in the chain in such a way that states are 
transferred with perfect (JJ-JH or with arbitrary high fidelity iflOl - lfTBI ; in addition, some minimal ex- 
ternal control on the chain dynamics was also introduced in order to achieve similar results |[T8l - l|26| . 
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Quite recently it was shown l27l that particular types of quasi-one dimensional uniformly coupled 
chains, can achieve perfect state transfer, provided that the natural dynamics of the chain be supple- 
mented with some global control pulses at regular intervals of time. Then it was shown in [27] and 
using a different geometry in [28], that perfect routing of states in higher dimensional networks, from 
any point to any other point, is also possible. The advantage of these schemes was that they allowed 
simultaneous routing of multiple states and also the possibility of overcoming some of the restrictions 
of the previous protocols, notably introducing some robustness to local imperfections in the network. 
It is important to note that the introduction of external control is not in contrast with the spirit of quan- 
tum state transfer through the natural dynamics of the chain, as long as the external control is global 
and does not address individual qubits in the network. 

While these schemes are to some extent robust against localized imperfections in couplings, i.e. 
by routing around known defects in the network, the new element of global control inherent in these 
kinds of schemes brings about the question of their robustness to inaccuracies and imperfections of 
external control. We can ask to what extent the fidelities of these schemes is affected by imprecision in 
the timings and the direction of the applied fields in the global pulses, which are necessary for perfect 
routing of states. 

It is on these novel stochastic disturbances that we shall focus in the present work. Specifically we 
focus on the quasi-one dimensional chain of fl27l . which is the basic element in higher dimensional 
uniformly coupled networks. We show that for a moderate value of the noise in global control, one 
can still achieve a high value of fidelity. Interestingly, we find that the transfer protocol appears to be 
less affected by correlated noise, as compared to un-correlated one; this result may be of great interest 
in the actual realization of realistic spin chain channels, since it suggests that externally induced time- 
correlations may protect the efficiency of the spin transmission lines. 

Although in general such noises affect the dynamics of the full chain in an analytically intractable 
way, we show that the fidelity of the state transfer can be exactly determined and analyzed. More 
specifically, we are able to compute the fidelity of the protocol, when the direction of the applied field 
and the timings of the pulses are not precisely tuned because of the presence of external stochastic 
noise. 

The structure of the paper is as follows. In the next section we give a brief account of perfect state 
transfer in the quasi-one dimensional chain of ll27l . Section 3 deals with the general formalism that 
allows to incorporate in the protocol the effects of the presence of disorder in the pulses. Section 4 is 
instead devoted to the actual computation of the fidelity both in presence of un-correlated disturbances, 
and in the more interesting case of correlated noise. The final section contains a brief discussion and 
outlook. 

2 Perfect state transfer in a uniformly coupled quasi-one dimen- 
sional chain 

The prototype of many-body system that has been used in many protocols is the XY spin chain, 
consisting of a linear array of N sites, to each of which a spin- 1/2 operator with Cartesian components 
X, Y, Z is attached. The dynamics is then described by a Hamiltonian of the form 




(1) 



2 



where the sum is over the various bonds in the array (see Fig.(l) for the specific example discussed in 
the following). This type of interaction preserves the total component of the spin along the z-direction 



0, (2) 



and moreover does not evolve the uniform background state where all the spin projections along the 
z direction are up, conventionally called the vacuum state |0): H\0) = 0. Further, let us denote with 
|i) the single excitation state corresponding to the situation in which all spins are up except the one 
at site i which is down; the collection { |i) }, i = 1,2, ... ,N, of these N states, form a basis in the 
single excitation subspace of the total Hilbert space of the system. 

In order to transfer an arbitrary qubit state ( Q ] from site 1 to site N, one may apply the following 



v/ 3 

simple protocol. First, embed the qubit into the initial lattice state a\0) + /3|1). Then, let it evolve 
according to the natural chain dynamics through the lattice until it becomes, at a suitable instant of 
time, the state a\0) + /3|N); from this state, one can then extract the original qubit from site N. 

As mentioned before, one way to achieve this perfect state transfer is to carefully engineer the 
coupling constants; this has been discussed in 0, where it was shown that a linear XY chain with 
local couplings of the form J„. rn = y/ n(N ~ n) 8 n +i,m can indeed perfectly transfer a qubit to the 
end of the chain, at the specific time t = tt. It is interesting to note that for N = 2 and N = 3, 
the couplings will be uniform and indeed it has been shown that these are the only uniformly coupled 
chains which can achieve perfect state transfer. 

The specific chain analyzed in ll27l is shown in Fig.(l). The chain is uniformly coupled in the 
sense that all the couplings have the same modulus. The presence of —1 couplings allows this chain 
to be broken up into a direct sum of sub-chains with just two and three sites, achieving perfect state 
transfer. To see this, consider the Hamiltonian pertaining to this chain; with reference to the labeling 
of the sites of Fig.(l), it can be expressed as 



H = (|1)<2| + |1)(3| + h.c) + (J2)(4| |3)(4| + h.c. ) + ( |4)(5| + |4)(6| + h.c. ) - 
which can be conveniently re-written as 



H = V2 



1}<(2,3)+| +h.c.) + |(2,3)_)(4| + \4)((5,6) + \+h.c. 



(4) 



where 



Ik) ± Ik - 

(k,k + 1)±) = ■ 



V2 

The orthogonality of the states | (k, k + 1 ) ± ) shows that the diamond-shaped lattice can be studied 
as a chain formed by two-site and three-site elementary components as shown in Fig.(l); we will refer 
to the latter as "virtual chain" and work with it from now on: it consists of K = (N — 4)/3 three-site 
sub-chains plus an initial and a final two-site sub-chain. A convenient basis, in the single-excitation 
Hilbert sub-space is then given by {\i)}, i = 1, 2, . . . , N, made of the following orthonormal vectors 

|1) = |1>, |2)= |(2,3) + ), |3) = |(2,3)_), |4) = |4),... (5) 

... |JV-1) = |(N-1,N)_), |JV) = |N). 

Note that the vector states |3i + 1), i = 0, 1, . . . , (N — l)/3, correspond to single spin down at the site 
3i + 1, like in the diamond-shape representation of the chain; however, unlike in that one, the other 
vectors represent a single-excitation being in a superposition of a spin down at one site and down at 
the subsequent one. From now on we shall work within the representation based on the orthogonal 
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Figure 1: (color online)Above: the quasi-one-dimensional chain introduced in [27|; solid lines indicate bonds 
with coupling equal to 1, dashed lines those with coupling equal to -1. Below: its equivalent description in terms 
of two- and three-site perfect state transfer chains. 



states non-boldfaced states in (0 corresponding to the disjoint union of small chains. 



As mentioned before, the time-evolution operator Ut = cxp(-itH) maps the single-excitation 
sub-space into itself. In the representation corresponding to the above orthonormal basis, it can be 
split into the following orthogonal sum of blocks 



tu?> 



V 



\ 



u 



(2) 



u 



(3) 



u 



(K+2) 



(6) 



The unitary blocks C/ t and {/ t , corresponding to the initial and final two-site sub-chains, can 
be represented as 



u 



(1) 



u 



(K+2) 



-itV2o 



cost-v/2 — isintv2 
— % sm 



(7) 
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by means of the Pauli matrix a x = ( ^ ^ ) , while the remaining ones can be represented as 



U 



-i2t 



1 



l + cos2t -iV2sin2t -l + cos2^ 



= - -i\/2sin2t 2cos2i -iy/2sin2t \ , j = 2, 3, . . . , K + 1 , (8) 



^-l + cos2< -i\/2sin2t l + cos2£ 
in terms of the x component of a spin-1 operator: 

1 /° 1 ° N 

S x = — 1 1 

v 2 \0 1 Oy 



(9) 



The use of a x and S x as generators of the blocks of the unitary time-evolution matrix Ut follows 
from the expression (01 for the Hamiltonian, whereby a x represents the first term in the sum, S x the 
second one and so on. One easily sees that, for t = n/(2y/2), respectively t = 7r/2, the unitaries 
in (O, respectively dH), realize prefect state transfers over the elementary two-, respectively three-site 
sub-chains. 

Besides the unitary evolution given by Ut, the transfer protocol introduced in 11271 involves in- 
stantaneous strong pulses P acting globally on the lower sites of the diamond chain in Fig.(l). The 
explicit form of the operator P is most simply given in the language of the virtual chain and in the 
basis (0: 

A \ 



P = 



1 



\ 



(10) 



V 



where a x appears K + 1 times, each one of them coupling \i) and \i + 1) for i = 3?i + 2 with 
n = 0, 1, . . . , K. Thus, except for the initial and final sites of the virtual chain, the pulse P transfers 
any single excitation at an end point of a given sub-chain to that of the contiguous one. The basic 
idea of the transfer protocol of |27l is to use the unitaries U t to transfer a single excitation along the 
sub-chains, while using P to make it jump from one sub-chain to the next. In practice, the composite 
dynamics of the system, starting from the initial time t = to the final time t = tt, is then described 
by the unitary operator 



~ U tf -t K+1 PU tK+1 -t K P ■ ■ ■ Ut 3 - tl PU tl 



(11) 



Suppose the times t l5 1 2 , ■ ■■, tx+i at which the pulses act are chosen as 

7T 7T 7T 



t 



1 2^2 



(12) 



Then, from the explicit form that ([6), (|7]i and ^ take in this case, one easily sees that the initial state 
a|0) + /3| 1), where |0) denotes the ground state, will be perfectly transformed into a\Q) + /3\N) at the 
final time 

7T 7T 
tf = — + -K . 

f V2 2 



(13) 



From this final state, the single qubit 1^1 can be extracted from site N, thus realizing its perfect 
transfer from the beginning to the end of the chain. 



3 Noise in the protocol 

As mentioned in the introduction, in order to transfer a generic qubit state 



a = cos - , 13 = sin -e~ # , 0<6><7r,0<(/><27r, 



(14) 



along the chain, one first embeds it into the left hand site of the chain itself as the state 

|tt)=a|O)+0|l) , 



(15) 
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where {|i)}j^ x are the basis vectors in (0 and |0) is again the ground state. Under the unitary action 
of the time-evolution (fTTT i. is transformed into \^tt) = Ut s |^); since H\0) = 0, 

\^ tf )=a\0)+p(j2Mtf)\^ ■ (16) 

The protocol purpose is to use the chain dynamics to transfer the state from site 1 to site N; the 
state of the qubit at site N is obtained by performing the partial trace (tr^) over single-excitation 
states involving all sites but the 7V-th one of the projection |^ r t/ )(^ r t/ |; this gives a 2 x 2 density 
matrix 

„/ fl fr \*( tMi6 ( t \\ f\a\ 2 + \f3\ 2 (l-\Mt f )\ 2 ) aF1> N (t f y\ 

p N {tf)-ir N Mtf))mi)\-[ a'PMtf) m 2 \^(tf)\y ■ ( } 

The robustness of the transmission along the chain can be measured by computing the fidelity of 
the final mixed state at the site N with respect to the initial pure one embedded at site 1; it is given by 

F V , = (V>|M*/)IV>) = M 2 + 2 \a\ 2 \(3\ 2 Re^ N (t f )) + |/3| 2 (|/3| 2 - \a\ 2 )\i> N (t f )\ 2 . (18) 

Note that it involves only the last coefficient ippf(tf) in the expansion (fl6i l. Its average over all input 
states, given by [T) 

F -=^J q d ^J Q ddsm9F 4 , = l + ±\ t jj N (t f )\ 2 + lRe^ N (t f )), (19) 

is a measure of the overall robustness of the transmission protocol. 

The fidelity is maximal, that is = 1, in the case of the ideal perfect protocol embodied by the 
unitary evolution lit # m CD with times as in ( TTZb and pulses as in ( TTOb . The success of the perfect 
state transfer depends on the precise control over the timing of the pulses and the direction of their 
field; in particular, the pulse should be exactly of the form ( TTOb to perfectly transfer the state from 
one sub-chain to the next. Also the timing of these pulses should be precise and synchronous with the 
timing required for perfect state transfer within the sub-chains. A disorder in either the field direction 
or timing of the pulses may drastically reduce the fidelity of the final state. 

In the following, we shall consider the case where external noise affects the protocol by modifying 
precisely these times and pulses. In practice, instead of those in (fT2l , we shall consider modified times 
of the form (as before, the initial time is set conventionally to zero): 

* 1 = 27i + T1 ' 

t j+1 = -^=+j J + T j+1 , 1.2 A. (20) 



2V2 2 

7T 7T 

V2 + 2' 



where r = {ri}^ 2 is a stochastic process with random variables Tj distributed according to joint 
probabilities Vtimeij) = Ftime(TK+2, Tfsr+i, ■ • ■ , Ti). The stochastic process may be stationary or 
not, correlated or not and the stochastic variables may take real values in a discrete or continuous state 
space. For sake of simplicity, we shall only suppose them to have zero mean-values. 

Analogously, we shall consider noisy pulses P(6i) that, while keeping the block form ( [Tol l, will 
no longer consist of local pulses represented by a x , but by 

P( 0. ):= ( isin ^ C0S M , ( 2i) 
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that reduces to a x in the limit 0, = 0; here, 9 = {^i}^ 1 is also a stochastic process with probability 
distribution V pu lse{6) = V pu ise(8K+i, Ok , • ■ • , 9\) consisting of real stochastic variables 9j with zero 
mean values. 

We shall collect the two stochastic processes into a single one, [i — (r, 9), with joint probabilities 
V(fJ.), that may even account for possible correlations between them. Then, in presence of such kind 
of noise, the unitary evolution ([TTT i will be replaced by 

= U(t k +2 - tk+i)P(9k+i)U{t k +i - tk) ■ ■ ■ U(r 2 - t 1 )P(9 1 )U(t 1 ) , (22) 

where the dependence on the stochastic variables in the various contributions is explicitly shown. 
Because of its very construction, despite the presence of a perturbing noise, U { t f will map the single- 
excitation sub-space into itself. Thus, for each realization of the noise, the initial state dT~5T > will be 
mapped into another single-excitation state 

TV 

\¥ t f) = u\f\*) =a|0) +/5 5>P 1) (*/)K> , (23) 

with a reduced density matrix at the YV-fh site, (tf) = tr fr\^^ (t f)) (t f)\, given again by 

the matrix in (fTTt . but with if>N(tf) replaced by ipjf 1 (tf), and similarly for the corresponding averaged 
fidelity (see (O). 

However, a physically meaningful state for the system can only be obtained by averaging over all 
realizations of the noise; the density matrix representing the state of the chain at the final time tf is 
therefore given by 

p{tf) = (p W it f )) = J dfiV(n) \®tf}{&tf\ , (24) 

where the integration region is given by the space spanned by the values that the stochastic variables 
(t, 6) can take. Then, the reduced TV-site state at time tf will become 

p N (t f ) = *rjv(p(t/)) = (trjv(p"(*/))) , (25) 
so that the fidelity averaged over the disorder and over all initial states will be given by 

(F) = J &»V{n)FM =\ + l(\^\tf)f) + l(Re(^\t f ))) . (26) 



4 Fidelity in the presence of noise 

We are now faced with the task of computing the fidelity (F) of the state d23l with respect to the initial 
state to be transferred. For this type of dynamics, it is obviously impossible to determine the explicit 
form of the final state given any arbitrary initial state. Nevertheless, as explicitly shown in (f26), in 
order to evaluate (F), only the coefficient ijjff (tf) in the expansion d23l is really needed. 

As already stressed, the crucial observation is that, in order to determine i^(tf), one needs 
just follow the change of the last non-zero entry of the vector |1) under the sequential action of 
operators of the form (O and ([H]). In doing so, we shall explicitly write only those components of 
the transformed vector affected by the various unitary blocks. Let us then consider the initial state 
|1) = (1,0,0...,0) T ; the action of 

, _ / -sin(\/2Yi) -i cos(v^n)\ 
{Tl) \-icosiV2n) -sin(V2n)J 
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transforms it into 



vector 




Then, the first pulse P(9\) maps the relevant two-component 
After that, the relevant three-component vector 



sin 8 1 cos(v / 2t"i) 
i cos^x cos(v / 2ti 



is turned by U (2 > (t2 — n) into 



/— i(l — cos(2(t 2 — ri))) cos 0i cos(v / 2t l ) n 
— I a/2 sin(2(r2 — ri)) cos 6>i cos(v / 2ri) 
z(l + cos(2(r 2 — 7"i))) cos 6*i cos(v / 2"r L ) 



(27) 



The last component corresponds to the basis vector |5); the second pulse P(02) turns it into the basis 
vector |6) multiplying it by cos 02- As such it is then subjected to U^ 3 \t3 — T2). Continuing in this 



way, the last non-zero entry of the final vector, that is the coefficient ip^ (t/) which we need, reads: 



(-1) 



K+l 



K+l 

J COS0j 



(V2( 



/K+l , 

'(v5n)(n- 

\ i=2 

c+2 - rRT+l)J 



+ COs(2(Ti - Tj-l)) 



(-1) 



Xw <Pn , 



(28) 



where x^v denotes the first bracket, namely the contribution from random pulses, while <p T N that from 
random time intervals between pulses. By averaging over the noise, one can then compute the mean 
fidelity (F) and thus study the impact of the noise on the robustness of the communication line. In the 
next subsections we will study in detail state transfer degradation along the chain both in presence of 
independent and correlated stochastic processes. As we shall see, the transfer protocol turns out to be 
more robust in the latter case, i.e. in presence of time-correlations. 



4.1 Independent noise 

The most common noise likely to affect spin chain communication lines is that generated by un- 
corrected disturbances: it can be described by totally independent stochastic variables {8i} and {rj} 
with uniform distributions. In this case, the probability density V(fJ,) = V(8, r) factorizes 

V{n)=V p ulse{e)Vtv m e{T) , (29) 
Ppulseid) = PpulseidK+l) Vpulse^K) ■ ■ - V pulse(®\) , (30) 
Vtimeij) = Vtime(jK+l) Vume{TK+l) ■ ■ •T'Wie^l) , (31) 

where, for simplicity, we have assumed the same probability distribution for all pulses and all time 
variables. Thus, the computation of (V'jy ) an d (I^tv^ 2 )' nee ded in the evaluation of the fidelity, 
simplifies, 

= {-l) K+l ( X 9 N ) m , (|^| 2 > = <(x^) 2 } (K) 2 > , 02) 

reducing d24t to the product of integrals for each stochastic variable. 

In addition, we shall assume the random variables to be uniformly distributed in the intervals 
9i G [—eg, eg] and n G [— e T ,e T ], around the perfect transfer values 6i — and r, = 0, so that 
V P uise(9i) = l/(2e#) and Vumeiji) = l/(2e T )- However, note that, besides on n, the quantity <j> T N 
in (l28l depends also on differences of the variables Tj. Therefore, in evaluating the averages (4> T N ) 

and ((t/)^) 2 ), it is convenient to introduce a new set of independent random variables, Si = Ti — Tj_i, 
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Figure 2: (color online)Average fidelity for various values of eg = e T — e in terms the chain length, (left 
figure, for a short chains of length up to ~ 60, right figure for long chains of length up to ~ 1200: N ~ 3K). 
From top to bottom, e — 0.01, 0.02, 0.04, and 0.06. The dashed line represent a fidelity of 0.95 and the dotted 
line represents the maximum fidelity, F = 2/3, reached using a classical communication line. 



i = 2, . . . , K + 2; being linear combinations of two uniformly distributed stochastic variables, these 
differences obey a triangular distribution 



2e T - \6j 
4e? 



-2e T <Si< 2e T . 



Taking this into account, the explicit computation yields: 

(\^\t f )f) = 



sin eg 


K+l 


sin(v / 2e T ) 


3 




e e 




\/2e T 




u 



sin 2e T 
2e T , 



K 




1 / sin2\/2e r 

2 1 1+ 2V2e T 



' 3 + 4 f^ 



sin4e T 
4e T 



(33) 



(34) 



K 

(35) 



All terms appearing in the square brackets above are less or equal to one, becoming smaller and 
smaller as eg and e T , measuring the strength of the noise, increase. Consequently, as the chain be- 
comes large, both (ip^\tf)} and (lip^itf)^} become small. Therefore, in presence of uniformly 
distributed noise, the fidelity (F) in (|26] i as function of the length K of the transmission line, decreases 
following a power law, the faster the more eg and e T differ from zero (see Fig. (2)): the advantage of 
using a quantum communication line with respect to a classical one becomes then rapidly ineffective. 
Nevertheless, if the errors induced by the noise are reasonably small (less than one percent), the aver- 
age fidelity remains above its corresponding classical value F = 2/3 for quite long chains (N > 900). 



4.2 Correlated noise 

The presence of correlations in the disorder affecting the spin chain is a concrete possibility in actual 
experimental realizations of the transmission protocol. Indeed, the presence of correlations among 
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subsequent pulses is likely to occur in practice due to the inevitable imperfections of the apparatus 
which produces them, as well as in the time intervals between pulse activations. In such instances, 
correlations are likely to arise, affecting in various ways the robustness of the spin chain transmission 
line. We shall concentrate on noise involving the pulses, i.e. the stochastic process 9; the disturbances 
affecting the time-intervals between pulses, described by the process r, can be similarly treated and 
will be briefly discussed at the end of the sectionQ 

In presence of correlations among pulses, the probability density P pu ise(9) can no longer be 
written as the product of independent probabilities as in d30l . The simplest generalization of <T30b 
involves joint probabilities with one-step correlations based on the conditional probabilities 

r pu ise(O l+1 \e i ) =q8(6 i+1 - 9i) + (1 - q)Vpuise(6i+i) , (36) 

specifying the probability that the stochastic variable corresponding to the i + 1-th pulse take the 
value conditioned on the stochastic variable corresponding to the preceding pulse having taken 
the value 6i, this being valid for alH = 1,2, ... ,K. The parameter q e [0, 1] measures the amount 
of correlations between the stochastic variables 6*; + i and 6*;, which is maximal for q = 1, while, for 
vanishing q, Qi+i and $i are independent stochastic variables. Then, we shall assume the stochas- 
tic process 8 = {9i} to be characterized by one-step correlations, namely that its joint probability 
distributions are of the form 

'Ppulse(0K+l, ■ ■ ■ ,02, Si) = VpulseiSK+l^x) ■ ■ ■ 'Ppulse (8 2 I #1 ) Vpulsei^l) ■ (37) 



In determining the fidelity in ( 1261 ). one now needs to use this expression in computing the averages 
over the noise. In practice, since the noise affects only pulses and not the time-intervals between them, 
the quantity d28l ) which enters d24l i together with its square modulus, reduces to (— 1) k+1 Xn> mus 
one needs evaluate integrals of the form 



d f(9K+l) V p ulse(0K+l\0K) ■ ■ ■ f{@2) Impulse N0i)/(0i)7V*e(0i), (38) 

where f{6) is either cos 9 or cos 2 9. 

In order to estimate the effects of correlated noise on the behavior of the averaged fidelity, we 
shall assume that the stochastic variables 9{ take only three possible values, —eg, and eg, with a 
probability distribution given by (0 < p < 1/2): 



V pu i se {9i) = { (39) 




Then, the Dirac delta in (1361 1 becomes a Kronecker delta and the integrals of the form ([38} reduce 
to sums: J d K+1 9 — » J^Ot e K +i an< ^ can ^ e cast ' n a com P act form by adopting a transition matrix 
formalism. That is, one introduces the orthonormal vectors 



eo) = , |0) = 1 , \eo) = , (40) 






a probability vector \Vpuise) with components (9\P pi ase) = r P pu ise(9), and collect the conditional 
probabilities into a 3 x 3 transition matrix with entries (9\¥ pu i se \9') = Vp U ise{d\9'). Explicitly, using 



4 Also here, for simplicity, the two processes 9 and r are assumed to be independent; correlations between pulses and time 
intervals are surely possible in principle, but certainly less likely than those between 0's and t's variables themselves. 
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( l36l and (l39l l. one finds: 




q + (i-q)p (i-q)p { 1 ~q)p 

| (l-g)(l-2jj) g+(l_g)(l_2p) (l-g)(l-2p) 

(l-g)p (i-g)p g + 

Further, by introducing the diagonal 3x3 matrix 

/(o) 

the average (l38l can be formally rewritten as the following matrix element: 

1= (e\(¥P pulse ) K ¥\v pu ise) , 

where the final vector (0| is the sum of the three basis vectors, explicitly given by 

(6|= Yl (°\ = ( 111 )- 

6 = {±e g ,0} 



(41) 



(42) 



F 



(43) 



(44) 



(45) 



By recalling that f(eg) is either cos eg or its square, this result allows evaluating for any fixed K 
the effects of correlated pulse noise in the fidelity d26l l as a function of the parameters q and p (see 
Fig. (3)). As expected, the efficiency of the qubit transfer through the chain degrades in presence of the 
correlated noise, but in a less dramatic way if compared with its behaviour in presence of disturbances 
with no correlations; indeed, almost perfect transfer is achieved for high correlated noise even when 
the error parameter p, the probability for the stochastic variables 6i to differ from the perfect transfer 
value 6i = 0, is as large as 1/2. 

In particular, for q = 1, i.e. when the correlations between subsequent pulses are maximal, one 
finds: 

(F) = 1 - p + | [2 (cos e e ) K+1 + (cos e e ) 2A '+ 2 ] , (46) 

clearly showing that external stochastic noise containing correlations, hence some sort of correlations, 
disturbs in a milder way the spin chain transmission protocol. Specifically, for long chains, as K 
becomes large, the averaged fidelity reaches the asymptotic value 



(F)~l-p. 



(47) 



which can still be close to unity, provided p is sufficiently small. 

This result has to be compared with the one obtained in the case q = 0, when no correlations are 
present and all stochastic variables 0i are independent. Also in this case, the pulse noise contribution 
to the averaged fidelity can be exactly computed, yielding: 



(F)= 1 - + 1 - 
W 2 6 



1 -4psin 2 (e e /2) 



K+l 



1 — 2psin 2 eg) 



K+l 



(48) 



Since < p < 1/2, the square brackets above are always < 1, so that as K increases the fidelity 
rapidly approaches its asymptotic value of 1/2; this is precisely the behaviour encountered in the 
previous section while discussing independent noise. 

Similar results are obtained when correlations are present in the stochastic variables r,;, affecting 
the time intervals between the pulses: the joint probabilities Vtimeiji\Tj) can be taken as in ( l36l ). 
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Figure 3: (color online) Average fidelity for various values of q in terms of the error probability p. From bottom 
to top, q = (no correlation), 0.9 and 1 (full correlation). For all the curves e — 0.5 and K = 100. 

By assuming that the stochastic variables Tj take only the three possible values — e T , and e T , with a 
probability distribution Vtime ( T i ) similar to the one in d39l l, the computation of the noise contributions 
to the averaged fidelity can be treated as in the previous case, leading to contributions of the form ( l44l ). 
However, note that now, except for the first contribution, the matrix F is no longer diagonal since it 
involves time differences; indeed, instead of F in (l43l l. one has to use one with entries: (t|F|t') = 
/(r — t'), where, recalling d28l ). f(c T ) is either (1 + cos 2e T )/2, or its square. 

The behavior of the averaged fidelity in terms of the probability p, for different values of the 
correlation parameter q, at fixed K and e T , is qualitatively similar to the one discussed before in the 
case of correlated pulses, given in Fig. (3). In particular, also in this case one observes that the fidelity 
is less affected by the presence of correlated noise, to the extent that when q = 1 it acquires a constant 
value, independent from the length of the chain: 



This result is easily understandable; indeed, when q = 1, all intermediate three-site sub-chains (lower 
picture in Fig.(l)) remain perfect state transfer chains even in presence of noise, and only the first 
and the last two-site sub-chain fail to transfer the state perfectly, so that the actual length of the 
transmission line becomes effectively irrelevant. 

5 Discussion 

We have studied the effect of imperfections in the external control in schemes for perfect transmission 
of quantum states through a quasi-dimensional chain ||27l . Such chains are to be connected to each 
other to form larger two and three dimensional networks 1271 , 1281 through which qubit states are to be 
routed from any point to any other point through the natural dynamics of the underlying XY Hamilto- 
nia when assisted by global control pulses. These protocols are by construction robust against known 




(49) 
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localized imperfections in the network structure (un-desired couplings, etc). Our study shows that as 
long as the quasi-one dimensional chains is concerned, such schemes are also robust against impre- 
cision in the sequence of applied global pulses, at least for moderate level of noise and for moderate 
lengths of the chain. Remarkably we have found that when the noise in successive applications of 
the pulses are correlated, the efficiency of the protocol is less damaged compared with the case when 
there is no correlations. Despite the complications of natural dynamics intervened by global external 
pulses, we have been able to derive exact expressions for the fidelity of state transfer, by taking ad- 
vantage of the sequential dynamics of the many-body state and following only the evolution of the 
relevant coefficient which is necessary for calculation of the fidelity. 
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